FIG. 9. Finite-size scaling curves of the resistivity, R{iijJn) I Rq, for the (2-|-l)-dimensional 
XY-model with random bonds correlated along the imaginary time axis with / = 1/2. From the 
top, the sets of curves represent the interpolated curves for K = 0.670,0.680,0.685,0.6875,0.690, 
and 0.695. We choose a = 0.40,0.35,0.30,0.27,0.24, and 0.20, respectively. The universal conduc- 
tivity determined from these curves is a* jaq = 0.49 ± 0.04. 
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FIG. 8. Finite-size scaling behavior of L-^Pxxi^) and (i^/i^)/>^^(0) for the (2+1)- dimen- 
sional XY-model with correlated random bonds along the imaginary time axis with frustration 
/ = 1/2 in the spatial plane. The transition point determined from these crossing curves is 
K* = 0.6875 ± 0.025. Inset: The scaling behavior of L^p^xi^) with respect to the single variable 
(K - K*)L^I''. K* = 0.6875 and 1/z/ = 1.0 are used. 
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FIG. 7. Finite-size scaling curves for the resistivity, R{iijJn) I Rq, for the (2-|-l)-dimensional 
XY-model with random bonds correlated along the imaginary time axis. From the top, the sets of 
curves represent the interpolated curves for K = 0.440,0.445,0.4475,0.450, and 0.455 respectively. 
We choose a = 0.88,0.55,0.39,0.28, and 0.16, respectively. The universal conductivity determined 
from these curves is a* jaq = 0.27 ± 0.04. 
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FIG. 6. Spin-spin correlation function of the (2-|-l)-dimensional XY-model along one spatial 
axis, say x, and along the imaginary time axis, (x) and ( + ) represent actual data points and 
the solid lines represent the interpolated curves. The parameters used in the interpolation are 
C'x = 0.2865, yx = 1.054, C-j- = 0.2975, and = 0.994. From the ratio of and y^, the dynamical 
exponent is determined to be 2; ~ 1.06. 
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FIG. 5. Finite-size scaling behavior of L-^Pxxi^) and (i^/i^)/>^^(0) for the (2-|-l)-dimensional 
XY-model with correlated random bonds along the imaginary time axis. The data points for 
the systems with non-integral size are obtained by the linear interpolation of the data obtained 
for systems which have nearby integer sizes. The transition point determined from these crossing 
curves is K* = 0.4465 ± 0.025. Inset: The scaling behavior of L^p^xi^) with respect to the single 
variable {K - K*)L^/''. K* = 0.4465 and 1/z/ = 1.0 are used. 
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FIG. 4. Finite-size scaling curves of a(iijJn)/ <^Q for the pure (2-|-l)-dimensional Xy-model. 
From the top, the sets of curves represent the interpolated curves for K = 0.704, 0.706, and 0.708 re- 
spectively obtained for different size systems. Here we choose a = 0.68, 0.61, and 0.48, respectively. 
From these curves, the d.c. conductivity is determined by the value at ujn/'^o — ct'^o/i^nLT) = 0. 
The universal conductivity obtained at the transition point is a* /gq = 0.52 ± 0.03 in this case. 
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FIG. 3. Spin-spin correlation function of the pure (2-|-l)-dimensional XY-model along one 
spatial axis, say x, and along the imaginary time axis for frustration / = 1/2. (x) and ( + ) represent 
actual data points and the solid lines represent the interpolated curves. Because of the frustration 
the spin-spin correlation function along a spatial axis oscillates slightly. The interpolated curve 
of Gxif') is the correlation function for spins separated by even numbers of lattice constants. The 
parameters used in the interpolation are Cx = 0.4495,^2; = 1.071, = 0.5341, and y^- = 1.075. 
From the ratio of and y^, the dynamical exponent is determined to be 2; ~ 1. 
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FIG. 2. Finite-size scaling behavior of L-^Pxxi^) and (i^/i^)/>^^(0) for the pure (2+1)- dimen- 
sional XY-model with / = 1/2. The transition point determined from these crossing curves is 
K* = 0.707 ± 0.001. Inset: The scaling behavior of L^p^xi^) with respect to the single variable 
(K - K*)L'^I''. Here we use K* = 0.7068 and 1/z/ = 1.5. 
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FIGURES 

FIG. 1. Probability distribution of energy in Monte Carlo simulation of the pure boson model 
for the case with / = 1/2 and / = 1/3 for different size systems indicated in the figure. The 
absence of an energy barrier separating two phases shows that the transition is continuous. 
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IjvK, 1.5. For the disordered case, we study the model for / = and / = 1/2. We obtain the 
quantum dynamical exponents z = 1.07ib0.03 and z = 1.14ib0.03, respectively. The critical 
exponent 1/;/ 1 is obtained, the value of which is consistent with Chayes lower bound 
condition [36]. The universal conductivity is (t j (Tq = 0.27 it 0.04 and (t j (Tq = 0.49 it 0.04 for 
/ = and / = 1/2, respectively. Disorder in the model has little effect on the value of the 
universal conductivity (at least at / = 1/2). Also it does not wash out the length scale set 
by the frustration, although the critical exponents are changed. 
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this model. Fig. 7 shows the curve ol the resistivity lor different couplings. Here we plot 
the resistivity instead ol the conductivity. In a disordered boson model, S0rensen et al. [23] 
could fit their resistivity data in a Drude form. Even though our data expressed in the form 
ol resistivity are not fit in a Drude form, better scalings are obtained Irom them than those 
obtained Irom the corresponding conductivity data. The curves are obtained Irom system 
sizes 10 X 10 X 9, 12 x 12 x 11, 14 x 14 x 13, and 16 x 16 x 15. Since the errors caused by 
the uncertainty ol the critical coupling are dominant compared with those caused by the 
uncertainty ol we use systems with nearby integer sizes. From those curves, we determine 
the universal conductivity ol the model: a* j gq = 0.27 it 0.04. 

With Irustration / = 1/2 in the spatial planes ol the model, we have z 1.14 it 0.03, 
K* = 0.6875 ± 0.025, ly = 1.0 ± 0.3 (see Fig. 8), and a*/aQ = 0.49 ± 0.04 (see Fig. 9). 
Again, the error bound on the values ol z and p are the rough range for good fitting. The 
correlation lunction calculation e 1.17. The universal conductivity in this case is also 

approximately doubled compared with the unlrustrated model. This means that the type 
ol disorder we consider does not wash out the magnetic length. As for the critical coupling, 
even though the shape ol the lattices used are slightly different in the simulations with 
and without randomness along the imaginary time axis, the critical coupling is consistently 
suppressed in the disordered case. This result is consistent with the recent simulation [29] 
ol a model high- Tc superconductor with columnar defects along the c-axis, where the Tc is 
increased because ol the vortex pinning effect due to the random disorder. 

IV. CONCLUSIONS 

We have calculated the universal conductivity for a random-// boson Hubbard model. For 
the luUy Irustrated (/ = 1/2) pure system, we obtain ct/ctq = 0.52ib0.03. The conductivity is 
almost doubled, as predicted by Granato and Kosterlitz [34], compared with the unlrustrated 
system. For / = 1/3, we find ct/ctq = 0.83 it 0.06 so that the conductivity is approximately 
tripled. The phase transition in the above two cases is continuous, and we find z ^ 1 and 
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tion of the correlation function, and we find Ijv 1.5. The universal conductivity in this 
case is approximately tripled compared to the unfrustrated case, which is consistent with 
the existence of three different ground state configurations. 

B. Random- [/ model 

For the (2+l)-dimensional XY-vaodeX with random correlated bonds in the imaginary 
time direction, we have tried several values for the quantum dynamical exponent, z. Since 
the lattice sizes along the imaginary time direction determined by Eq. (11) are in general non- 
integer numbers, we adjust the aspect ratio to make the lattice sizes as close to integers as 
possible. More than 1000 disorder configurations are taken near the transition. When / = 0, 
we have the best scaling behavior 1.07. Fig. 5 shows that, in this case, the scaling 

curves of LrPxx{0) and [L^ / Lr)pTT{0) for systems with the sizes of 10 X 10 X 9.070, 12 X 12 X 
11.023, 14x 14x 13,andl6 x 16 x 14.997 cross around a transition point K* = 0.4465 ±0.025. 
Since the stiffness depends weakly on the aspect ratio, these curves were obtained by linear 
interpolation of the data from the systems which have sizes corresponding to the nearby 
integers. Here the aspect ratio is As = 0.772. For significantly different either we cannot 
obtain the crossing behavior at all or L^PxtX^) and [L'^ / L-r)pTT{^) cross at different points. 
The error bound for the number z^ which gives good scaling, is roughly z = 1.07 it 0.03. A 
slightly larger number was obtained by Lee et al. [29] in a similar model at / = 1/4, however 
their error bound is too big to meaningfully compare the numbers. The value of z is also 
consistent with the direct calculation of correlation functions as shown in Fig. 6, where we 
find z 1.06. The value of ;/ = 1.0 it 0.3 is obtained from the scaling curve shown in the 
inset of Fig. 5, where the error bound is the rough range for good scalings. In Fig. 5, we 
used Ijv = 1.0 and K* = 0.4465. This number for v is consistent with the simple extension 
of the Chayes lower bound condition [36] for the model with correlated disorder, which tells 
us that p > 2/d where d = 2 is the number of dimensions in which the disorder occurs. The 
error bound on p obtained in this fitting, however, is too big to confirm the condition for 
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an energy barrier between the two wells separates two phases. In finite systems, the energy 
barrier might be disguised by a finite size effect that will be reduced as the size of system 
increases so that the energy barrier increases as the size of system increases. We do not have 
any energy barrier increasing as the size of the system increases, and the result supports a 
continuous phase transition for both / = 1/2 and / = 1/3. 

Fig. 2 shows the scaling behavior of L^PxtX^) and [L'^ / L-r)pTT{^) for the pure (2+1)- 
dimensional XK-model with frustration / = 1/2 for systems of different sizes. We assume 
that 2; = 1, which is supported by the correlation function calculation (see Fig. 3). The 
curves cross around K* = 0.707 it 0.001. The slopes of the curves give us information about 
the correlation length critical exponent p. The inset of Fig. 2 is the finite-size scaling of 
LtPxx{0) with respect to the single variable [K — K*)L^^'^. With l/f = 1.5 it 0.3, good 
scaling is obtained. This value of p is very close to p in the unfrustrated three-dimensional 
XY model [33]. The scaling behavior of the conductivity is shown in Fig. 4. From these 
curves, we obtain a* j (Tq = 0.52 it 0.03. The error range is mostly due to the uncertainty in 
the critical coupling, K* . Compared to the unfrustrated case where a* j gq = 0.285 it 0.02 
[12], the universal conductivity for the case / = 1/2 is almost doubled. This result was 
predicted by Granato and Kosterlitz [34] based on a picture of coupled XK-models [35]. For 
the fully frustrated case (/ = 1/2), there are two degenerate ground state configurations. 
Thus when we construct an effective Ginzburg-Landau action, we need two complex fields 
(or two species of bosons) fiuctuating around each ground state configuration. If we neglect 
the couplings between fields, we have the Gaussian model. The universal conductivity is 
proportional to the number of boson species in the Gaussian model. Since the 1/A^-expansion 
calculation of the universal conductivity [12] supports the notion that the Gaussian model 
is good as a first approximation, the universal conductivity of the fully frustrated case will 
be estimated to be doubled compared to the unfrustrated case. Our calculation supports 
this conjecture. 

We also study the case with / = 1/3, obtaining that K* = 0.841 it 0.03 and 
(T* I(Tq = 0.83 lb 0.06. Again z is assumed to be 1, which is consistent with a calcula- 
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where p^xi^^ L-^ / L^) is an universal number, and A is a non-universal constant. Similar 
behavior holds for [L'^ / L-r)pTT{^)- We use these properties to find the transition point, the 
quantum dynamical exponent, 2;, and the critical exponent, in the next section. 
The conductivity is obtained through the Kubo formula [12] 

(T^^ = 27raQ hm (21) 

where ctq = (2e)^//i. It is a (scaling) dimensionless quantity and has a universal value at 
the transition. In practice, however, the finite size contribution, which gives the residual 
superfiuid density, must be subtracted to find the universal behavior of the conductivity 
at the transition. Because of the frequency dependence and the finite-size contribution, we 
expect the conductivity h scaling form [12] 

^ ' = c{—-a — 7-) + ... (22) 

where loq is a given frequency such as the ultraviolet cutoff frequency. We set loq = 2tt j a in 
this work, where a is the lattice constant. We take the parameter a which gives the smallest 
deviation among the interpolated scaling curves for the different size systems. 



III. RESULTS 

A. Pure model with / = 1/2 and / = 1/3 

First, we consider the case of the pure boson model with frustration. We first have to 
check if the transition is first order or continuous in this case, because the scaling argument 
which supports the universal conductivity is based on the assumption of the continuous 
transition. A recent Monte Carlo calculation of the case / = 1/6 on a triangular lattice 
suggests that the transition is first order [30]. We use the Lee and Kosterlitz method [31] to 
check the nature of the transition at larger values of /. Fig. 1 is the probability distribution 
of energy, p(e), for / = 1/2 and / = 1/3, where p{e) = e~^^^/Tr{£;}.e~^^^ and e is the energy 
per volume. If the transition is first order, we expect a double-well type curve, in which 
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G,{r) = C,{r-y^ + {L,-r)-y^) (15) 

where Cx and Cr are fitting constants. Tlie results are discussed in the next section. Also 
one might expect the magnetic field perpendicular to the two-dimensional system to change 
the exponent z. We find that the external magnetic field has little efi'ect on z for the pure 
model and only weakly changes z in the random-// model. 
The frequency-dependent stiffness is given by 

(16) 

with 

= ^cos(AX-A,^) (17) 
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Miu^) = J2sin{AJr - A^^)e'-"^ (18) 

T 

where the frequency Un is the wavenumber in r-direction. In Eq. (16), (. . .) represents the 
thermal average and [. . .]avg represents the ensemble average over the impurity configura- 
tions. In this work, the contribution of the last term in Eq. (16) drops out because Lf is 
an integer in the pure model or because there is vortex-antivortex symmetry [i.e. f = or 
/ = 1/2) in the disordered model. 

The finite-size scaling ansatz for Pxx{^) is given by [12] 

PxAO) = ^P.AtL'^'',Lr/L^) . (19) 

where t = K — K* [K* is the critical coupling), p is the correlation length critical exponent 
((^ ~ t"''), and pxx is a universal scaling function. Near the transition, 

LrP..{0) = P..{0, LJL') + AtL'/'' + ... (20) 



In order to determine the different sizes of lattice witli tlie same aspect ratio, tlierefore, 
information about tlie value of z is prerequisite. For tlie pure boson Hubbard model with 
the particle-hole symmetry, it is known that z = \ [15]. Quenched disorder or long-range 
interactions, which alter the diverging behavior of correlation lengths, will change z [15,18]. 
For the model we are considering, we expect that impurities hinder the correlation along the 
spatial directions so that the correlation length along the spatial directions will diverge more 
slowly than that along the imaginary time direction. This certainly implies that z > \. The 
microscopic particle-hole symmetry property of the model (djidQj —djidO^ in Fq. (6)) 
implies that the occupation of the lowest delocalized eigenstate [i.e. the mobility edge state) 
does not change the total number of particles at the brink of the superconductor-insulator 
transition. Since the compressibility h form 

K ~ r' , (12) 

this condition implies that z < d where d is the spatial dimension. Thus the model is 
different from that with random on-site chemical potential with a short-range interaction, 
where it is expected the compressibility of the system is finite even at the transition and 
z = d. For the model we consider, therefore, we expect 

Kz <2 . (13) 

Due to the lack of the information about an exact value of z for the model we are studying, 
we have tried several numbers for z. We also find z by directly measuring the correla- 
tion functions. The correlation functions along one spatial direction, say x, and along the 
imaginary time direction have asymptotic behavior [10,32] 

Gr{r) ~ r-y^ (14) 

at the transition point, respectively, where y^^ = d — 2 -\- z -\- rj and = {d — 2 -\- z -\- r])/z. 
The value of z is obtained from z = y^Jy^- However in finite size systems with periodic 
boundary conditions, instead of Fq. (14) we use the form [32] 
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II. PROCEDURE 



The phase transition of the model defined in Eq. (8), is studied by Monte Carlo sim- 
ulations. We take simple cubic (2-|-l)-dimensional lattices whose sizes are systematically 
changed to fit data using finite-size scaling. Periodic boundary conditions are imposed and 
the ordinary Landau gauge is used for the phase change due to the magnetic field. We 
change the coupling strength parameter K to tune the transition of the model. 

The stiffness to a twist at the boundaries is measured to find the transition point and the 
critical exponents. The stiffness along a spatial direction is proportional to the superfiuid 
density, and has a finite value in the superfiuid phase but vanishes in the insulating phase. 
In a finite size system, however, it vanishes smoothly and there is a residual superfiuid 
density at the transition point, whose magnitude depends on the size of the system. From 
the size dependency, therefore, we can extrapolate to find the asymptotic scaling behavior 
by changing the system size while keep the same geometry of the system. Since the model 
is anisotropic in the presence of the correlated disorder along the imaginary time dimension 
or the frustration in the spatial planes, the diverging behaviors of correlation lengths near 
critical point are not isotropic in general, but they are related through the dynamical scaling 
hypothesis 



where z is the quantum dynamical exponent. This implies that in order to simulate the 
model to find the asymptotic scaling in a finite size system in which correlation lengths are 
limited by the size of the system, we need to take lattices with 



where As is the aspect ratio, and L and L-j- are the sizes in the spatial directions and in the 
imaginary time direction, respectively. In this work, we use lattices with = Ly = L and 
L-r which is determined by Eq. (If). 



(10) 
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k ^ y^e'"^' (5) 
if no >> 1. The model then reduces to the quantum rotor model 

n = J2 -P'f-T.J cos(^, - e, - a,) (6) 

where J = ^/n^to- Here we have shifted rii — no — > n^. When // = 0, we have the correspond- 
ing (2-|-l)-dimensional classical action, S[6]^ defined by 

Z = Tr^eye-^^ = Tr^eye-'^'^ , (7) 
which is obtained by the path integral transformation [28]: 

S[e] = - J2[Kr{i) cos(A^^,) + K cos(A^^,) + K cos(Aj,^,)] (8) 

r 

where i and r denote position on the spatial plane and the (imaginary) temporal axis 
respectively, r G (z, r) is a point in the (2-|-l)-dimensional lattice, and A-rOr = Or+i- — ^r, etc. 
Kr{i) is the coupling parameter along the (imaginary) temporal axis, depending only on 
the spatial position. Here we restrict the model to the particle-hole symmetric case // = 0, 
since otherwise the action is complex. The randomness of Ui in Eq. (6) is transformed into 
the randomness of K^- The frustration /, of the model is defined by the summation of the 
phase change around each plaquette so that 

27r/ = + Ajk + Am + Au . (9) 

This model is equivalent to a model high- Tc superconductor with columnar defects along 
the magnetic field recently studied by Lee, Stroud, and Girvin [29]. 

In Section H, we describe the properties of the model, including scaling hypotheses 
necessary to analyze the numerical data. We then study the pure model, in which Kr{i) = K , 
with frustration f=l/2 and f=l/3 in Section HI. A random-// model where are random 
numbers bounded by < j K < 2, with frustration f=0 and f=l/2 is also studied. Section 
IV is reserved for the conclusions. 
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In this work, we study the transition in a "phase-only" modeL Here we argue that 
fluctuations of the magnitude of the superconducting order parameter are irrelevant near 
the transition because they cost more energy than the fluctuations of the phase. This picture 
implies that there are coarse-grained boson islands in which lots of bosons condense to have 
well-deflned phase angles locally. These islands are connected through inter-island hopping. 
In this picture the phase angle is the only dynamical variable and the number operator and 
the phase operator are conjugate variables. The transition is determined by the competition 
of the density fluctuations and the phase angle fluctuations. Strong disorder and interactions 
in general suppress the density fluctuations locally and cause more fluctuations in phase 
angle. The external magnetic held causes a frustration in the phase couplings in this model. 

More explicitly, we consider a lattice boson Hubbard model with random on-site inter- 
action: 

^ = E - ^0 - fif - + t,^b]k) , (1) 

where bi (b]) is the boson annihilation (creation) operator at site z, rii is the boson number 
operator at site z, Ui is the interaction between bosons at site z, tij is the nearest neighbor 
hopping matrix element, and no-\- fi (no is an integer and —1/2 < // < 1/2) is the parameter 
which determines the average background boson density. In the presence of an external 
magnetic held, we have 

= toe'^'^ (2) 

27r f^i -> 
= — / A{x) ■ dx 

where $o is the flux quantum and A[x) represents the vector potential 
well deflned phase angle 6i at each site z, 

- iA 

~ I do, 

and 
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(3) 

When there is a 
(4) 



amplitudes [17]. In particular, the conductivity at the zero-temperature superconductor- 
insulator transition has been predicted [18] to be a universal number, which is analogous to 
the universal jump of the superfluid density at Kosterlitz-Thouless transition. Experimen- 
tally measured values of resistivity (or resistance per square) at very low temperature in 
the vicinity of the superconductor-insulator transition are close to the quantum resistance 
[1,4-6], Rq = /i/(2e)^, in various systems, even though no consensus concerning whether the 
temperature is low enough to reach the critical regime and whether the measured values are 
truly universal [i.e. independent of material parameters) has emerged yet. Theoretically a 
model of interacting bosons moving in the presence of disorder has been proposed to capture 
the appropriate universality class [15,18]. Considerable theoretical effort has been devoted 
to the study of the phase transition in this model by Monte Carlo simulations [19,20] and by 
renormalization group calculations [21,22]. Also the universal conductivity was calculated in 
the various versions of the model through different techniques [12,23-27]. Even though those 
calculations support the existence of a finite critical conductivity at the transition, the cal- 
culated values are either inconsistent among different studies or a factor of 2-3 away from 
experimental values. So far, theoretical studies mostly concentrated on the case without 
magnetic field. Here we want to study the phase transition and the universal conductivity 
in an interacting boson model in the presence of an external magnetic field. We also consider 
the effect of disorder given in the form of a short range repulsive random interaction. 

This zero-temperature phase transition can be understood by adopting the charge- vortex 
duality picture [11]. Near the transition, vortices and anti-vortices are induced by quantum 
fiuctuations. In the superconducting state, delocalized bosons condense whereas vortices 
are localized, yielding no dissipation. On the other hand, in the insulating phase, bosons 
are localized and vortices condense. According to this scenario, at the superconductor- 
insulator transition both bosons and vortices are expected to be mobile [12], yielding a finite 
conductivity because the motion of charged bosons carries current while that of vortices 
produces voltage. The external magnetic field changes the density of vortices and breaks 
time reversal symmetry of the vortex distribution thereby modifying the universality class. 

4 



I. INTRODUCTION 



Through recent experimental studies in various two dimensional systems [1-9], the pres- 
ence of a zero-temperature phase transition between the superconducting phase and the 
insulating phase has been convincingly suggested [10]. The transition has been tuned by 
changing parameters, such as the thickness of homogeneous films [1-3], the strength of 
charging energy in Josephson junction arrays [4], the strength of disorder [5], or the mag- 
nitude of the external magnetic field applied perpendicular to the two dimensional systems 
[6-9]. We view this transition as a boson localization problem, in which electron pairs are 
treated as point bosons small on the scale of the diverging phase correlation length near 
the transition. This picture is certainly suitable in granular superconductors and Josephson 
junction arrays where the electron-pair wavefunctions are well defined in each grain and the 
superconductivity is obtained by the establishment of a long-range phase correlation among 
them. Recently it has been proposed [11,12] that even in homogeneous films, the zero- 
temperature superconductor-insulator transition will be controlled by phase fiuctuations of 
the electron-pair wavefunction. The picture is based on the assumption that below a mean- 
field temperature, electrons form pairs but the true superconducting temperature is lowered 
by the enhanced phase fiuctuations in low dimensional systems. The size of the electron- 
pairs is finite at the transition. Hebard and Paalanen [7,13] reported some evidence that 
in homogeneous films, such as composite InOx films, the magnitude of the superconducting 
order parameter is finite at the transition, supporting the idea that the relevant fiuctuations 
at the transition are the long wavelength twists in the phase of the pair wavefunction. We 
cannot rule out the possibility, however, that in these systems the transition might possibly 
be induced by the vanishing of the magnitude of the order parameter [14]. 

Based on the assumption that the transition is a continuous phase transition and that the 
phase fiuctuations of the superconducting order parameter give the primary contribution to 
the singular part of the free energy, scaling theories have predicted some universal properties 
of the transition such as critical exponents [15,16] and dimensionless combinations of critical 
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the case / = 1/2. 
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Abstract 

The universal conductivity at the zero-temperature superconductor-insulator 
transition of the two-dimensional boson Hubbard model is studied for cases 
both with and without magnetic field by Monte Carlo simulations of the 
(2-|-l)-dimensional classical XY-model with disorder represented by random 
bonds correlated along the imaginary time dimension. The effect of magnetic 
field is characterized by the frustration /. From the scaling behavior of the 
stiffness, we determine the quantum dynamical exponent z, the correlation 
length exponent z/, and the universal conductivity a*. For the disorder-free 
model with / = f/2, we obtain 2; ~ f, f/z/ ~ f.5, and a* /gq = 0.52 ± 0.03 
where gq is the quantum conductance. We also study the case with / = f/3, 
in which we find g* /gq = 0.83 ± 0.06. The value of g* is consistent with 
a theoretical estimate based on the Gaussian model. For the model with 
random interactions, we find z = f .07 ± 0.03, z/ ~ f , and g* /gq = 0.27 ± 0.04 
for the case / = 0, and ^ = f .f4 ± 0.03, z/ ^ f , and g*Igq = 0.49 ± 0.04 for 
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